Vibrational coherence in electron spin resonance in nanoscale oscillators 
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We study a scheme for electrical detection, using electron spin resonance, of coherent vibrations in 
a molecular single electron level trapped near a conduction channel. Both equilibrium spin-currents 
and non-equilibrium spin- and charge currents are investigated. Inelastic side-band anti-resonances 
corresponding to the vibrational modes appear in the electron spin resonance spectrum. 
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Inelastic effects arising due to coupling between 
charge carriers and local vibrational modes (vibrons) in 
nanoscale electronics devices has gained an enormous in- 
terest recently. Peaks and dips often observed in the 
differential conductance of molecular electronics deviceJ^l 
may indicate strong effects from electron-vibron cou- 
pling. Effects from vibrons have been investigated^ ^ in 
molecular quantum dots and single electron transistor, 
in Josephson junctions,^ and on surfaces using scanning 
tunneling microscopy (STM).^ 

The interplay between vibrons and charge carrier is 
expected to generate dynamical signatures also in the 
spin current or spin-dependent transport. Such dynam- 
ics should consequently be observable in electron-spin 
resonance (ESR), which will thus allow for electrical de- 
tection of both spin an d vi bron modes. In this paper, 
we apply the ESR set-ufP^ to a molecular quantum dot 
with electron levels coupled to vibrons, and we show 
the emergence of anti-resonances in the spin-current at 
frequencies equal to integral numbers of the vibrational 
mode. The anti-resonances can be explained as inter- 
ference between opposite spin tunneling electron wave 
functions traversing different molecular excitations. Such 
information would be useful not only to the conventional 
semi-conductor industry, but also to novel research direc- 
tions such as spintronics and molecular electronics. Dis- 
tinct from earlier work on ESR set-up we are address- 
ing the novel signatures arising from the electron-vibron 
coupling, which will manifest in both equilibrium spin 
current and nonequilibrium spin dependent transport. 

We model the resonator to be oscillating with fre- 
quency luq, where the vibrational motion is weakly cou- 
pled to the electrons with strength A. We consider the 
dynamics of a single molecular level Eq coupled to exter- 
nal thermal baths. The level is spin split by the exter- 
nal magnetic field Bq, — = uj^ = q^bBq, where 
g and are the gyromagnetic ratio and Bohr magne- 
ton, respectively. The spins are coupled by a rotating 
magnetic field i3i(coswit,sina;if) applied perpendicular 
to i?o, and we assume 2gfiBBi <C LUo,ujr- We employ the 
model H — Hc + Ti-d + "^t, where 

T-Ld = ^[e^ + A(a^+a) + J7nff/2]n<, 

-g^iBBi{d\die''^^^ + dU^e-''^'') + ujoa^a, (1) 



describes the molecular states, while Tic = X^fecr ^knka- 
and Ti.T — X^fecr ^kc\,jda + H.c, are the Hamiltonians for 
the bath and the tunneling, respectively. Here and 
dj. create an electron with spin cr =t,i in the bath and 
molecule, respectively, Uka- = c^crCfeo- and no- = d^^dcr. 
The operators and a denote creation and destruction 
of the vibrational mode. 

We transform the system into the rotating reference 
frame of the magnetic field through Hrf = e^''fTle~^^f + 
i{dte^''f )e~^''f , with the unitary transformation Srf = 
—i{u!it/2)[ni — + X]fe('^'=i ~ "-fct)]' ™ order to elimi- 
nate the time-dependence from the Hamiltonian at the 
cost of introducing a shift in the electronic energies, i.e. 
e)^^ — efe + crwi/2 and e]j^ — + CT(jJi/2, where the factor 
CT = ±1. The spin- split of the conduction channel elec- 
tron energies originates from the magnetic pumping field 
through the hybridization between the localized level and 
the conduction band. The pumping propagates energy 
from the molecule to the conduction channel and gen- 
erate the spin chemical potentials — —oLOxjl (with 
reference to ep — 0) in the conduction channel. The 
frequency of the oscillating magnetic field can thus be 
regarded as the (spin) bias applied to the system. De- 
spite the spin-imbalance, however, the charge chemical 
potential is still ^ = (/Z| -I- /ij^)/2 = 0.1^ 

Although the system itself is to be considered in equi- 
librium, the one photon imbalance between the spin 
channels generates a non-equilibrium condition for the 
two spin projections of the electrons. An electron in the 
spin down channel can thus tunnel into the local spin 
down level e^i^ . The rotating magnetic field flips the spin 
projection of the localized electron and thereby the elec- 
tron can tunnel into the spin up channel, and a stationary 
current builds up by repeated tunneling. 

The coupling between the vibrational and electronic 
degrees of freedom is de-coupled by the canonical trans- 
formation = e^p'^Tlrfe~^p''' with Sph = (A/wo)(«^ ~ 
a) Ua-. Through this transformation the energy levels 
of the localized states are turned into eg. — e^J — \^ /luq, 
while the charging energy U = U — 2A^/wo; and the tun- 
nelling Hamiltonian is changed into ^ka-('^kiTc\.^d„X + 
H.c), where X — exp [— (A/tJo)(a^ " a)]. In the present 
study we assume a weak coupling between the electronic 
and vibrational degrees of freedom, and that the spin- 
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currents through the system are small. It is then justified 
to neglect narrowing eS'ects on the tunneling between the 
conduction channel and the molecular levelP 

In the atomic limit and U — 0, the molecule is reduced 
to a simple driven two-level system. It is characterized 
by a coherent weight transfer, Rabi oscillations, between 
the two spin states, which is complete at resonant rotat- 
ing frequency toi = uj^. The spin oscillation period of 
T = 27r/f], where n = A'^ + 4{gfiBBiy is the Rabi 
frequency and A = uji — Ur denotes the detuning from 
the resonance. We transform the molecular electron op- 
erators by 



u = 



cos (p — sm ( 
sin d) cos 6 



(2) 



where tan0 = 2g^BBi/{il — A). The molecular elec- 
tronic states are diagonal in the new representation giv- 
ing the molecular Hamiltonian E^rcl^Ccr + tjc^^c-^c'i^ci, 
with ^ (et + q - (Tn)/2 ^ Sq ~ X^/ujq - cr^l/2. 
The spin- (7 current I„ is preferably written as 

= ^tr y r,{/<,HG>M + U{Lo)]G<{Lo)}dLo, 

(3) 

where r^- = u^^r, T = 2T:Y,k \vk\'^5{uj - eu) and 



„ „ I sin^ (h sin 6 cos ( 

Ut = T^Uit'^, Ui = • ^ J, 2 a. 

>- ' >- \ sm cos ffl cos 6 



(4) 



where is the y-component of the Pauli matrix, whereas 
fa{^) = .f{^ ^ fJ-a) is the Fermi function for the spin a 
channel. 

The current contains the lesser (greater) Green 
functions (GFs) G<(» = {0^^^}^^'- They can 
be calculated using G<(» = G''S<(»G", where 
e.g. the retarded GF is defined through G^^,{t) = 
{—i)dit){{ccr{t),cl^,{0)}) and similarly for the advanced 
one. The canonical de-coupling procedure of the electron- 
vibron coupling casts the GF into the product of an elec- 
tronic and vibronic part as 

= G:At){X{t)X\0)),,,. (5) 

with da{t) = e*'"='*c^e-^'^'='*, X{t) = e*^-i.*Xe-*^""*. 
The renormalization factor caused by the electron- vibron 
couphng is calculated as {X{t)X'^{0))yib — e~*^*\ where 



$(t) = (A/c^o)2[nB(l-e^"«*) + (nB + l)(l- 



-iujQt 



)], with 



riB — {e^^° — 1) -'^.'i^We then calculate the electronic 
GF G^^, in the mean field approximation, in which the 
Kondo resonance effect is neglected. For arbitrary on- 
site charging energy, its Fourier transform is given by 
G^„g{uj) — and GJJ.g. — G"^, where 



iJ - - {1 - {n^))U 



{lu-E„ + iT /2){u; -E„-U)+ i{ns)UT /2 



(6) 



and {n„) = Im j G<{uj)duj / {2n). We then find 

n 

(7) 

where In{z) is the nth modified Bessel function and z = 

2[\/ujQfyJnB[nB + l\- 

In the case of weak electron-vibron coupling the con- 
tributions to the self-energy S from the electron-vibron 
interaction is negligible, hence, the lesser (greater) self- 
energy can be approximated by 

i:< = if^T^ + if^T^, s> = -zrt-zri + s< . (8) 

This leads to the spin-cr current — eF^ J T[/o-(ijj) — 
fff{oj)]dLLj/h, where the transmission coefficient 

T^\G^uj)~Gl{uj)\^sm^(l)Cos'^(t). (9) 

We notice that the transmission T is equal for the two 
spin channels, which is expected in the stationary regime, 
thus giving the total spin-current 1^ — '''aa-^'y ~ 2^t- 
The form of the transmission coefficient suggests that 
the spin-current can be interpreted as an interference be- 
tween tunneling electron wavefunctions in the conduction 
channel, which are coupled by the molecular level. 

The interference interpretation is especially appealing 
in the context of electron-vibron coupling. For simplicity 
consider the case of vanishing effective charging energy, 
U — {U ^ 2X^/ujq), although the arguments remain 
true for arbitrary tj. Then, the first factor in Eq. ^ 
can be written as (setting cj^^" — E^ ^iT /2) 

r2/„(z)e"^"«/2 



\Gl 



E 



[ljJ — ojI — nujQ)(u! — id'', — ncuo) 



(10) 

The main ESR peak is given at uji = ujr such that 
fl = 2gfiBBi and E^r = Eq — X^/loq — gq^bBi, corre- 
sponding to the n = term in the transmission coeffi- 
cient. Because of the electron-vibron coupling, additional 
features in the spin-current are expected to occur at fre- 
quencies LJi = ujr + nujQ, corresponding to the vibrational 
side-bands. Due to the destructive interference between 
tunneling electron waves passing through different con- 
duction (spin) channels, these satellites to the main ESR 
peak appear as dips in the spin-current rather than peaks. 

In order to illustrate this argument, we consider the 
first two terms in Eq. ( 10 1, that is the terms with n = 0, 1 



n'[I^iz)\G\iLo)Gl{Lu)\' + ltiz)\G\{u; - u^) 



xGliu; - tJo)Pe^"° + 2ReIo{z)h{z)G^uj)Gl{uj) 



xG?(w - Wo)G?(t^ - ujo)e^'^<'/'^] 



(11) 



where the first two terms add positively to the transmis- 
sion, and peak at w = E^r and uj = Ect + luq, respectively. 
The last term, proportional to 



Re 



-Re 



n/{ujo + iV) 



{uj- E^ -Luo- iT/2){uj - Ei+ iT/2) 

n/jmp + iT) 

{lu- E^-ujQ- iT/2){uj - E^ + iT/2) 



(12) 
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is negligible at uji = cUr since then E-^ 
to that the two contributions cancel each other 
+ ujq, on the other hand, we have _B| + ojq s; 



El , which leads 



As uji 
E\ , since 



~ Wo- Therefore, 
roughly equals 



the first contribution in Eq. (12) 



i/[i + (r/c.o)^] 

{uo-E,Y + {Tl2y 



(13) 



while the second contribution is negligible. The expres- 



sion in Eq. (13 1 peaks around lo — E\^ w ii'i + wq, and 



contributes destructively to the total transmission coeffi- 
cient in Eq. (11). An estimate of the ratios between the 



third and first, and third and second terms in Eq. 
at Wi ~ uJr + i-^Oi yields the lower bounds 



(11) 



h{z) 



2/o(z) 



2h{z) 



-'3-»/2£(c^o), (14) 



respectively, where C{uJo) = + (wo/r)^] = r2/[l -f 

(r/wo)^]- We, thus, find that the transmission is sig- 
nificantly reduced when the detuning A equals the first 
vibrational side-band. Including the remaining contribu- 
tions to the transmission, i.e. summing over all n, provide 
similar reductions in the transmission at all frequencies 

iOi= IjOr + nujQ. 

We calculate the spin current by solving Eq. ([6| self- 
consistently. The equilibrium spin current through the 
molecular level is plotted in Fig. [T] as function of the ro- 
tating frequency wi, illustrating the main ESR peak at 
uji = LOr and the vibrational anti-resonances at cui = 
ujr + nujQ, n {ujo/ujr = 1/2 in the plot). At van- 
ishing correlation energy, the spin-current decreases for 
increasing coupling strength A, which is understood as 
an effect of the density being distributed among an in- 
creasing number of vibrational side-bands for increasing 
electron-vibron coupling, c.f. Fig. [2]ja). Increasing spin- 
current for increasing correlation energy can be explained 
by the same effect, see Fig. ^h). 

It is easy to generalize the above theory to two leads 
and non-equilibrium conditions. The voltage between the 
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Vdj =0.35 U,.'M =100 






leads is eV = hl — IJ-r, and in each lead we have the 
spin-imbalance such that = (/j,^| -I- ^^|)/2, x = L,R. 
The current /^o- for the spin-cr current flowing from the 
left lead into the molecule is written as (with obvious 
notation) 



ha = 



h 



J r^{r^T,[/i,(c^) - fn,{uj)]+T'^%[,fLa{^) 



- hsioj)] + T^'TsihA^) - fBAuj)]}doj, (15) 

where % = |G|Cos^0-|- C^j^ s\n^ , whereas 7^ is the 
transmission coefficient given in Eq. Here, also 

f^a{Lo) — f{uj — /ixtr)- The expression for the current 
in Eq. ( 15 1 is obtained by the observation that the lesser 
(greater) self-energy in this case is given by 



(16) 



We identify the contributions in Eq. ( 15 1 by the first 
being the usual charge transport as derived by Meir and 
Wingreen,^^ the second contribution is the one discussed 
above in Eq. (|3| , and the third contribution accounts for 
the spin current between the leads. 

The charge current between the leads, Ic = J2cr^L(j, 
becomes 

ic = ^Y.J r^r«(r, + %)[fu^) - /^.(u;)]^^^, (i7) 

which is just the sum of the different transmission con- 
tributions between the leads. As one would expect, 
Ic lacks the interference effects that occur in the spin 
current, which becomes clear by noticing that 7^ + 
Ts = IG^pcos^^-l- \Gl\'^ sin^ (j). The spin-current Is = 

J2a '^ajLa IS givCU by 

Is = I yr^{2r^T,[/it(c.)-/ii(c.)] + r^r,[/it(c.) 



fR]{t^) + fRli^^) - /Li(w)]}dw, 



(18) 



which contains three contributions. The first contribu- 
tion has the same origin as discussed above in the sin- 
gle medium case; the second contribution accounts for 




FIG. 1: (Color online) Equilibrium spin-current dependence 
of A and U. Here, cUr = 2, gfisBi = 0.2315, T = 4^/5/25, and 
fcflT — 10, in units of ljq. 



FIG. 2: (Color online) Local molecular DOS for X/loq = 0.35, 
and U — (a) and U /luq = 10 (b). Other parameters are as 
in Fig. [1] 
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FIG. 3: (Color online) Bias voltage dependent spin-current 
for A/wo = 0.35 and (7 = 0. Other parameters as in Fig. [l] 



the spin current between the leads; the third contribu- 
tion stems from the spin imbalance in the charge current, 
which arises from the spin-biased leads. 

The non-equilibrium spin current is plotted in Fig. [3] 
showing its dependence on the bias voltage. For low ro- 
tating frequency, the spin current is dominated by trans- 
port that is assisted by the rotating magnetic field, c.f. 
first and second terms in Eq. (18 1, which provides the 



main ESR peak and vibrational anti-resonances analo- 
gous to the equilibrium case. Increasing frequency ui 
increases the potential barrier for a molecular level spin- 
flip. Hence, ac magnetic fleld assisted transport becomes 
suppressed, in analogy with the equilibrium situation. 
The non-equilibrium conditions do, however, enhance 
tunneling between the leads of electrons that does not 
undergo spin-flips when in the molecule, i.e. the contri- 
bution from the third term in Eq. 



18 1 increases. 



We have, for simplicity, neglected effects on the vibra- 
tional coherence from the environment, which is justi- 
fled whenever the dwell time of the localized electrons 
= r ^ ciJo.^^ By studying the vibron mode life-time 



Tv, to the second order in the electron- vibron coupling A, 
we find ~ X'^T/[tt{luq + F^)]. In the present study 
we thus have ojqTv ^ 10^, for typical electron-coupling 
strengths and couplings between the local electron with 
the conduction channel, which justifies the introduced 
approximations. Assuming a vibron mode ~ 1 /^eV, 
provides a vibron life-time at least in the order of tens of 
ns, which should be sufficient for measurements. 

The herein reported anti-resonances are expected to 
occur more generally within the ESR set-up. In fact, 
we have studied the occurrence of the anti-resonances 
in systems where the local level is coupled to a general 
system with two or more levels, and we find that the 
anti-resonances will occur whenever the subsystems have 
direct interactions with one another. The nature of the 
interactions may be e.g. tunneling. Coulomb, spin-spin 
exchange interactions between electrons, or, as discussed 
in the present paper, fermion-boson interactions. 

In summary, we have studied a scheme for electrical 
detection, using ESR, of vibrational coherence in molecu- 
lar single electron level trapped near conduction channel. 
We have shown that the electron-vibron coupling gener- 
ates anti-resonances in the spin current at frequencies 
equal to integral numbers of the vibrational mode. The 
anti-resonances can be explained as interference between 
opposite spin tunneling electron wave functions travers- 
ing different molecular excitations. Observations of the 
vibrational anti-resonances does not require extremely 
low temperatures, and since we using realistic parame- 
ters for in study, we believe that our findings should be 
within experimental reach. 
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